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ABSTRACT 

The geometric calibration of the Planck satellite is investigated, specifically those 
parameters which require the use of the science data for their extraction. Methods for 
the recovery of these geometric-calibration parameters from point source detections in 
the science data are presented, together with the accuracies which may be achieved. 
These methods apply to the a posteriori evaluation of these parameters using all the 
mission data, and may be incorporated into the initial stages of the construction of 
the Planck final compact source catalogue. It is found that this process achieves the 
pointing requirements, provided that the geometric-calibration parameters vary only 
slowly in time. Indeed the errors in the pointing reconstruction due to the geometric 
calibration parameters may be seen to approach those due to the star tracker. 

Key words: cosmic microwave background — techniques: miscellaneous — cata- 
logues — submillimetre — radio continuum: general — galaxies: general 



1 INTRODUCTION 

Planck is a European Space Agency satellite designed to 
produce high-resolution temperature and polarisation maps 
of the CMB. It is equipped with detectors sensitive to a 
wide range of frequencies from 30 to 857 GHz, split be- 
tween two instruments, HFI and LFI, the high and low 
frequency instruments, respectively. The full-sky coverage 
and frequency range provided by Planck will also enable 
the construction of a unique compact source catalogue. Be- 
low 100 GH z it will be a si g nifica nt improvement over that 
of WMAP, feennett et alJ 12003ft - ). in terms of sensitivity 
and hence in the number of sources detected, while above 
100 GHz it will be the only full-sky survey for many years 
to come. 

This paper investigates the possiblity of recovering the 
geometric-calibration parameters, defined in section |3~T1 as 
part of the initial stages of the construction of t his source 
catalogue. Previous work, llHarrison et alJ ll2004ft h on the 
reconstruction of the geometric-calibration parameters con- 
centrated on the recovery of the boresight and roll angle 
parameters during the course of the mission; the recovery of 
the reference phase was not discussed, nor has it been dis- 
cussed elsewhere. While the monitoring and evaluation of 
these parameters is crucial over the lifetime of the mission, 
it is expected that much higher accuracies are achievable 
with an a posteriori evaluation of these parameters using the 
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entirety of the mission data. Hence it is expected that the 
construction of the final compact source catalogue, FCSC, 
will provide the definitive values of the geometric-calibration 
parameters. 

An overview of the FCSC is presented in Section |5] 
The geometric-calibration parameters are defined in Sec- 
tional together with a discussion on their use in the point- 
ing reconstruction of the Planck field-of-view and the ac- 
curacy requirements of the pointing reconstruction are as- 
sessed. The simulations generated to assess the performance 
of the geometric-calibration methods, outlined in Section ^] 
are discussed in Section The results of this analysis are 
presented in Section |SJ where we show that it is possible 
to recover these parameters to the required accuracies us- 
ing solely the bright extragalactic point sources, provided 
there is only a slow linear variation in the value of these 
parameters over the course of the mission. 



2 AN OVERVIEW OF THE FINAL COMPACT 
SOURCE CATALOGUE 

The Final Compact Source Catalogue, FCSC, should not be 
confused with the Early Release Compact Source Catalogue, 
ERCSC. The ERCSC will be released approximately 1 year 
after the start of the routine observations. Its primary pur- 
pose being to allow rapid follow-up observations from Her- 
schel and ground-based millimetre instruments. This cata- 
logue will be produced using the frequency-channel maps 
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from the data corresponding to the first full-sky coverage, 
this is expected to be the first ~ 7 months worth of data. 
The requirement on the release date, however, curtails the 
focus of the ERCSC to the brightest sources. Further details 
of the ERCSC implementation plan may be found in lGangal 
(2002). The FCSC will comprise of the full mission data for 
Planck and be a significant improvement over the ERCSC. 

The primary goal in the construction of the FCSC is 
to ensure the internal consistency and accuracy of the cat- 
alogue and final frequency maps. This goal has motivated 
the division of the catalogue construction into four major 
stages. 

(i) Stage fa: The accumulation of the detections from 
bright point sources in the ring data, where the ring data 
corresponds to a single spin axis positioning, is used to find 
positions for these point sources and any deviations from 
the current bestfit values for the geometric-calibration pa- 
rameters. It is this stage which forms the main focus of this 
paper. 

(ii) Stage fb: Solar system objects with their highly accu- 
rate positional data also allow the recovery of the geometric- 
calibration parameters. Given the relative motion of these 
bodies and the Planck satellite the methods employed are 
some what different from those in Stage fa. Hence the dis- 
cussion of the methods used in Stage I b is reserved for pa- 
per II, The Geometric Calibration of the Planck satellite us- 
ing Solar System Objects 

(iii) Stage 2 & 3: These stages involve the treatment of 
the intermediate-source detections on rings and the source 
detections in maps, respectively. It is envisaged that for the 
source detections from maps, it may be necessary to return 
to the ring data to ensure that maximum accuracy and con- 
sistency will be attained. The discussion of the methods in- 
volved in these stages of construction is left to paper III, on 
the final compact source catalogue 



3 GEOMETRIC CALIBRATION 

Planck will be inserted into a Lissajous orbit around the 
second Lagrange point of the Earth-Sun system, spinning 
about its axis once per minute. The line of sight being al- 
most perpendicular to the spin axis, hence the detectors 
will almost follow a great circle. The spin axis will nomi- 
nally be repositioned every hour, and the roughly 60 or so 
circles corresponding to a single spin axis positioning may 
be binned together to form a ring. The nominal spin axis 
passes through the centre of the solar panels and is directed 
away from the sun, thus maintaining the rest of the satellite 
in a cone of shadow produced by the solar panels. The scan- 
ning strategy is determined by the sequence of the nominal 
positions of the spin axis over the course of the mission. 

The geometric calibration is the process whereby all the 
line-of-sight positions of the detectors are recovered, at any 
time during the observations. The relationship between the 
pointing and ti me may depend on multipl e parameters, dis- 
cussed in full in Ivan Leeuwen et al.l i2002t) : the star tracker 
will provide many of these parameters. However, as will be 
discussed in Section 13.21 there is an uncertainty between 
the line-of-sight of the star tracker and that of the focal 
plane, this will produce offsets in the values of a few of the 
geometric-calibration parameters. These offsets may only be 




Figure 1. The figure shows the two angles which describe the 
position of the line-of-sight of a detector, LOSd, with respect to 
the spin axis position, SA, and the North Ecliptic Pole, NEP. The 
angle between the spin axis and the LOSd i s given by the opening 
angle to the detector, ad- The opening angle to the detector from 
the spin axis position describes the path followed along the sky by 
the LOSd- The position of the LOSd along this path is given by 
the phase, ipd- The zero point for the phase may be determined 
by the intersection point of the great circle which connects the 
spin axis position and the NEP and the path of the scan circle. 
This great circle may be defined as the reference meridian. 

recovered with the science data, and it is those geometric- 
calibration parameters which require this calibration with 
the science data that we discuss here. 

3.1 The Geometric-Calibration Parameters 

The position of the line-of-sight of a detector with respect 
to the spin axis position may be described in terms of two 
angles, which are shown in Figure □ The first, the opening 
angle, is the angle between the spin axis position and the 
line-of-sight of the detector. The opening angle defines the 
path described by the detector, the scan circle. The second 
angle, the phase, defines the position along the scan circle 
from a given reference point. This reference point may be 
given by the intersection point of the scan circle and the 
great circle connecting the spin axis position with the North 
Ecliptic Pole, NEP. This great circle, shown in Figure^ will 
be referred to as the reference meridian. Each repositioning 
of the spin axis defines its own reference meridian. 

If we assume the focal-plane geometry, as defined by 
the relative locations of the detectors to the position of a 
fiducial reference point, FRP, defined somewhere central in 
the focal plane, remains fixed, then instead of requiring two 
angles per detector we may describe the positions of all the 
detectors in terms of just three angles. These three angles, 
shown in Figure are the opening angle to the FRP, other- 
wise known as the boresight angle, ot frp , the phase ipFRP 
of the FRP, determined from the reference phase, lp re f, as 
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Figure 2. The figure shows the geometric-calibration parameters 
which allow a description of the position of the field of view, 
FOV, with respect to the spin axis position, SA. The boresight 
angle, appp is the angle between the FRP and the spin axis. The 
rotation of the focal plane around the FRP with respect to the 
nominal scan direction is given by the roll angle, p. The reference 
phase, ip re f, is the value of the initial phase at the point at which 
the FRP crosses the reference point, as defined by the intersection 
point of the reference circle and reference meridian closest to the 
NEP. The initial phase is measured from the first point observed 
on the reference circle. 



measured from the reference point defined by the reference 
meridian, and the roll angle, p, which is given by the rota- 
tion of the focal plane around the FRP, relative to a nominal 
scan direction defined within the focal-plane geometry, as il- 
lustrated in Figure El 

The position of the line-of-sight of a detector may now 
be described in terms of the nominal values for these angles, 
the phase of the d th detector is given by: 

Ipd — TpFRP + Xd(p) (1) 

and similarly the opening angle to the detector by: 

ctd = cxfrp + yd{p) (2) 

where Xd and yd are respectively, the scan and cross-scan 
positions of the d th detector with respect to the FRP and 
may be determined by the focal-plane layout and the roll 
angle. The scan and cross-scan positions of a detector for a 
given roll angle, p, may be found using the position given 
for the detector in the focal plane, (xd,o,yd,a)> which corre- 
sponds to the scan and cross-scan positions of the detector 
in the case of zero roll angle: 

( Xd \ — ( cos P sin P)( x d,o\ ,o\ 
\Vd J ~ y-sinp cospj\yd,oj l > 

The actual values of the geometric-calibration parame- 
ters, however, may deviate from their nominal values, as will 
be discussed in Section f3. 21 It is therefore useful to express 
the offset in the phase, Aipd, and opening angle, Aad, to a 
detector in terms the differences between the nominal and 
actual geometric-calibration parameters using equations 
toU 
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Figure 3. The figure shows the positions of the line-of-sight of the 
Planck detectors, filled circles, relative to the FRP which is shown 
by the cross. The black circles are the positions for the nominal 
scan direction, whereas the grey circles illustrate the effects of 
a non-zero roll angle on the detector positions. The dotted lines 
enclose detectors belonging to the same frequency channel. 



Aip d = 5ip F RP + x d (p) (cos(<5p) - 1) + j/d(p) sin(<5p), 
Aa d = SapRp - Xd(p) sin(<5p) +ya(p) (cos(r5p) - 1) . (4) 

The dependence of the recovered position of the line-of- 
sight of a detector may now clearly be seen in equation 0] 
the position along the scan is affected by the reference phase 
and the roll angle, whereas the position in the cross-scan di- 
rection is affected by boresight and roll angles. Information 
on the scan-phase correction Aipd may be obtained from the 
measurements of ipd for point-source transits. The correction 
to the opening angle Aad may be derived from the distribu- 
tion of intensities of the point-source transits. These correc- 
tions may be used to evaluate the offsets in the geometric- 
calibration parameters, this analysis is the subject of Sec- 
tional 

3.2 Pointing Reconstruction 

The pointing reconstruction of the satellite is achieved by 
the recovery of the the relationship between the line-of-sight 
of the star tracker and that of the focal plane. 

It is helpful to define two reference frames for Planck, 
the Satellite Reference System (SRS) and the Inertial Refer- 
ence System (IRS). The origin of both systems is the centre 
of gravity, CoG, of the satellite. The SRS is aligned with 
the principle axes of the satellite, and is shown in Figure [I] 
The x-axis is defined by the nominal spin axis which passes 
through the centre of the solar panels, and the CoG, and 
is directed away from the Sun. The z-axis is orthogonal to 
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Figure 4. The figure shows the Satellite Reference System, SRS, 
which is defined in terms of the principle axes of the Planck satel- 
lite. The x-axis is aligned with the nominal spin axis, the line-of- 
sight of the FRP lies in the z,z-plane, and the y-axis completes 
the right-handed triad. 



the x-axis and lies in the plane denned by the i-axis and 
the projected line-of-sight of the FRP. The y-axis then com- 
pletes this right-handed triad. It is necessary to define a sec- 
ond reference system, the Inertial Reference System, IRS, in 
which the inertia tensor is diagonal. The x-axis is now the 
actual spin axis, which is defined by the largest principle 
axis as determined by the principle moments of inertia of 
the satellite. 

These two reference systems are connected by a time 
varing rotation matrix which must be recovered in flight. As 
the mission progresses and consumables are depleted this 
affects the inertia tensor of the satellite and hence the axis 
of rotation, changing the relationship between the SRS and 
the IRS. The star tracker should recover this time varing 
relationship. However, the exact relationship between the 
star tracker reference system and the SRS is not known. 
This lack of an exact relation between the star tracker and 
the SRS and IRS will produce uncertainties in the recon- 
struction of the actual line-of-sights for the detectors. These 
uncertaint i es hav e been estimated to be ~ 1.3 arc-minutes, 
IChlewickH i2004) . and may be described in terms of offsets 
in the geometric-calibration parameters, which must be re- 
covered from the science data in order to meet the accuracy 
requirements discussed in Section f3. 31 

The relationships between the focal plane and the star 
tracker to the SRS, although nominally constant may indeed 
have a slow variation with time due to thermal effects. The 
offsets in the geometric-calibration parameters which must 
be established from the science data, may therefore be time 
dependent. Additional time variation may arise if the star 
tracker fails to completely establish the time dependence 
between the SRS and the IRS. The expected magnitude of 



these effects are unknown, but are thought to be less than 
1 arcmin. 

A more detailed discussion of the re ference systems and 
the at titude analysis may be found in Ivan Leeuwen et alJ 




3.3 Pointing Accuracy Requirements 

Figure|S|was created in order to assess the accuracy to which 
the pointing must be recovered to avoid compromising the 
reconstruction of the Ct power spectrum, which describes 
the magnitude of the Auctions in the CMB as a function of 
angular scale and is the primary science goal of the Planck 
satellite. It shows the effect of unknown random pointing 
reconstruction errors on the recovered Ci values. The solid 
grey curves correspond to the error on individual multipoles, 
whereas the dashed grey curves correspond to the error on 
multipole bins of 50. The shape of the grey curve is deter- 
mined at low multipoles, on the left of the figure, by the 
sample and cosmic variance. On the right at high multi- 
poles, however, the error is due to noise and the finite size 
of the beam. The decrease in the sensitivity of the beam to 
the higher multipoles may be accounted for at the cost of 
the exponential increase in the noise at the high £ values; 
as the noise is undiminished by the beam size. The series of 
black curves in Figure |S] show the effect of unknown random 
pointing reconstruction errors. The pointing uncertainties 
result in an effective smearing of the beam, resulting in a 
larger effective beam. This produces the reduction in the 
sensitivity to the higher multipoles, as seen in Figure [5] 
A more detaile d discussion of how th is figure is obtained, 
may be found in lHarrison et all i2004) . Ideally, the errors in 
the reconstruction of the pointing should be such that the 
additional errors introduced are less than that of the un- 
subtractable noise. This places an upper limit on the total 
pointing error of 9 arcseconds. 

It should be noted that correlated errors in the pointing 
will not produce effects as large as those for pointing errors of 
the same magnitude which are random. The effect of smear- 
ing the beam predominately in one direction is to increase 
the sensitivity to higher t multipoles as compared to the 
uniform smearing shown in Figure The pointing accuracy 
requirement generated by assuming unknown random point- 
ing errors should therefore exceed any requirements found 
from assuming correlated errors, and will represent the tight- 
est constraints on the pointing accuracy required. 



4 METHODS FOR THE GEOMETRIC 
CALIBRATION 

As discussed in Section [3.21 there will be an offset between 
the line-of-sight, LOS, as determined by the Star Tracker 
and the actual LOS. It is this offset which produces the 
systematic offsets between the actual and nominal values of 
the reference phase, boresight and roll angles, which may 
only be recovered using the science data. 

The methods for evaluating the geometric-calibration 
parameters therefore, need only assess the offset in the value 
of the parameters from their nominal values. These methods 
should also be able to cope with a slow variation in the 
values of these offsets over the course of the mission, due 
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Figure 5. The figure shows the unsubtractable errors on the 
Ces which are dominated by cosmic variance to the left and the 
finite beam size and noise to the right. The solid grey curves 
enclose the region corresponding to the error on each individual 
multipole, whereas the dashed grey curves correspond to the error 
on multipole bins of width 50. The black curves correspond to the 
additional errors on the reconstruction of the Cgs due to unknown 
random pointing errors. A beam of FWHM 5', with a noise of 
10.3/xK per beam and a fractional sky coverage of 0.7 have been 
assumed. 

to thermal effects, as discussed in Section f3.2l These offsets 
must be recovered to accuracies which will allow the total 
pointing error to remain below 9 arcseconds, as discussed in 
Section 13.31 

As discussed in Section loTTl offsets in the reference phase 
and the boresight angle produce errors in the recovered po- 
sition which are orthogonal, and hence may be solved for in- 
dependently. The offset in the roll angle, however, produces 
errors in position, in both the scan and cross-scan directions. 
In practice, however, the roll angle may be solved together 
with the reference phase independently from the boresight 
angle, as will be discussed in Section f4. II The calibration of 
the opening angles to the detectors and hence the boresight 
angle is discussed in Section f4.2l 

4.1 Reference Phase and Roll Angle 

The position of a point source in the scan direction may 
be constrained by the phase of the detection. Detections on 
multiple non-parallel scans which correspond to the same 
point source, may therefore be used to constrain the posi- 
tion of the point source on the sky. This process is illustrated 
in Figure [(J where the solid arrows represent the scan direc- 
tions. The position of the point source may be found from 
the intersection point of lines orthogonal to the scan direc- 
tion extending from the phase of each detection. Figure |S| 
shows these lines for the cases of no offsets in phase, dashed 
lines, and unreconstructed time dependent offsets, dotted 
lines. As seen in Figure|S|any unaccounted for offset in phase 
will affect the recovered position of the point source. Offsets 
in the opening angles to the detectors, however, will not 
affect the recovered position. Offsets in the opening angles 
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Figure 6. This figure shows the evaluation of the position of a 
point source, with, dashed lines, and without, dotted lines, in- 
cluding corrections to offsets in the reference phase or roll angle. 
The minimisation of the residuals between the expected phase of 
the source given its position and the observed phase of the de- 
tection allows an assessment of the offsets in the reference phase 
and roll angle as a function of time. 



only affect the angular separation between the recovered po- 
sition and the apparent position of the scan. The offsets in 
the phase of the detections, therefore, allow an offset in the 
roll angle to be discovered independently of an offset in the 
boresight angle. An offset in the roll angle is equivalent to 
a different reference phase offset for every detector, hence 
an offset in the reference phase may be distinguished from 
an offset in the roll angle. The minimisation of the resid- 
uals between the phase observed for the detection and the 
phase expected for the point source based on its recovered 
position, therefore provides a mechanism whereby the off- 
set in the reference phase and roll angle may be assessed. 
The methods presented here are similar to those emp loyed 
in the sphe re reconstruction for th e Hipparcos satellite. lESAl 
J 19971) and lLindeeren et alJ lll992h . 

The residual phase, Atpij, for the i th detection corre- 
sponding to the j th point source, may now be defined as the 
difference between the expected phase for i th detection given 
the position of the j th point source, ipij, and the observed 
phase of the detection, ipi. 

Aipij = tpij - ipi (5) 
The expected phase, ipij, may be evaluated using: 

Ipij = i>FRP i:j + Xd(pi) (6) 

where Xd is the scan phase of the detector, in which the i th 
detection occurs, with respect to the FRP and is dependent 
on the current value of the roll angle at the time of this 
detection, pi. The expected phase of the FRP for the i th 
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detection of the j point source, ipFRPi - , may be calculated 
using: 



tan(-0_f 



sin(AsA i — Xj) 



cos 0sa z tan/3j - svci0 S a z cos(X S a z - Xj) 



(7) 



where (Asa*, PsAi) is the spin axis position corresponding 
to the i th detection and (Xj,/3j) is the position of the j* 
point source. An initial position for the j* point source, 
(Xj,j3j), may be found by minimimizing the residual phases 
for all the detections of the j point source, hence solving: 



^ =0 and ^ A r 



dXi 



90, 



= 



(8) 



where the summation is over all the detections of the j th 
point source and Stpi is the residual in the phase of the i th 
detection which may be found using the relationship be- 
tween small changes in the phase of a detection and the 
resultant change in position, (8Xi,8f3i) : 

( sin PsAi — cos Qi sin Qi 



sx, 



( cosffgAj sin(A S A i - Xj) 
\ sin 2 at 



50, 



(9) 



where a.i is the current value, at the time of the i th detec- 
tion, of the opening angle to the detector in which the i th 
detection occurs, and SXi — Xj — Xi and 8/3, = 0j— it where 
(Xi,0i) is the position of the detection. 

Once an initial position for the j th point source has been 
found, the expected phase of the detection, ipij, and hence 
the residual phase, Atpij, may be evaluated. The residual 
phase, however, may also be related to changes in the ex- 
pected phase due to changes in the position of the j th point 
source, and the residuals in the phase of the detection, Aipd'- 

Atpij = Stpij + A^ d - (10) 

Using equation 0] and under the assumption that the offset 
in the roll angle is small, Atpd may be expressed as: 



Aip d = —5ip re f(ti) + yd{p) Sp(ti) 



(11) 



where the offsets in the reference phase and roll angle, may 
be expressed as a constant offset and a term defining the 
linear drift in time, 



Slp re f(ti) = (l/Ve/o + VWi*»)> 

8p(ti) = (p + p\U) . 



(12) 



Equation I1UI may now be written in terms of changes 
to the point source position and the offsets in the reference 
phase and roll angle, using equations^ and equation llll 

Alpij = AijSXj — Bij60j - (iprefo + ifrnfiU) 

+Vd(p) (po + piU) 

where, 

/ sin 0s Ai — cos on sin 0j ' 
'i = ( * ■ 2 I 



Bi 



cos0 S Ai sm(X S Ai - Xj) 



(13) 



A similar equation may be written for every detection, giv- 
ing i equations and 2 x j + 4 unknowns. The corrections to 
the point source positions and the values for the systematic 



offset and drift in the reference phase and roll angle may 
then be extracted by a nonlinear least squares analysis. 



4.2 The Boresight Angle 

As discussed in Section f3.il the opening angle to a detector 
is given by the angle between the line-of-sight of the detec- 
tor and the spin axis position. The opening angles to the 
detectors are determined by the focal-plane layout and the 
boresight and roll angles, hence if the offset in the roll angle 
has been accounted for, and the focal-plane layout is known, 
any remaining offset in the opening angles will be due to an 
offset in the boresight angle. 

Once the values of the reference phase, roll angle and 
the point source positions have been evaluated, as above, the 
positions of the point sources may be used to evaluate the 
value of the ordinate for each detection, where the ordinate 
is angular separation, in the cross-scan direction, between 
the point source and the path of LOS of the detector. 

The detections corresponding to the k th transit of the 
focal plane, in the cross-scan direction, of the j point 
source are accumulated for each of the detectors. It is then 
possible to find the value of the ordinates which correspond 
to the peak of the transits, in the cross-scan direction, for 
each detector in which the k th transit was observed. If there 
is no offset in the opening angle to the d th detector the value 
of the ordinate at the peak of the transit will be zero. Hence 
the offsets in the opening angles may be determined, at the 
epoch of each transit, by the analysis of the cross-scan tran- 
sits of the point sources. 

The offsets in the opening angle to the detectors for each 
transit of the focal plane may then be used to find the offset 
in the boresight angle at the epoch of the transit, 5aFRp{tk)- 
As with the reference phase and roll angle parameters, this 
offset may be expressed as a constant offset and a linear drift 
in time: 



5aFRp(tk) = 5aFRP + SaFRP^k 



(14) 



where tk is the epoch of the ft* transit. By assessing the 
value of SaFpp at every available epoch, epochs at which the 
focal plane transited a sufficiently bright point source, the 
time variation of the boresight angle may be investigated, 
and hence the systematic offset, SaFRp , and linear drift, 
SciFRPi, may be evaluated. 



4.3 Evaluating the Geometric-Calibration 
Parameters 

The evaluation of the reference phase, roll angle and point 
source positions is an iterative process, with offsets found 
from their inital values used in the further refinement of 
these parameters until convergence is reached. The offsets in 
the boresight angle should only need to be evaluated once, 
as the point source positions may be found independently 
of any offset in the boresight angle, as discussed above in 
Section 14.11 In practice, however, any significant offsets in 
the geometric-calibration parameters will affect which de- 
tections are classed as belonging to a single point source. It 
is therefore necessary to iterate over all the parameters, and 
reassign the detections as improvements in the values of the 
geometric-calibration parameters occur. 
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5 SIMULATIONS 

In order to assess the performance of the methods developed 
here in the reconstruction of the geometric-calibration pa- 
rameters, it is necessary to simulate a list of point source de- 
tections for the Planck mission. This requires an input point 
source catalogue, together with a scanning strategy and in- 
formation on the beams of the Planck detectors. Through- 
out this paper we have made the simplifying assumption of 
Gaussian beams. 



5.1 The scanning strategy 

While it is anticipated that this method will be applicable 
to any scanning strategy in which the circles corresponding 
to a single spin axis position may be binned together as a 
ring, only two potential scanning strategies for Planck were 
investigated here, a sinusoidal and a precessional scanning 
strategy. Where the sinusoidal scanning strategy may be de- 
scribed by: 



Afc = Ao + kd 
(3 k = Asm(n a X k ) 

and the precessional scanning strategy by: 
v = (Ao + k8) 



(15) 



sin(/3 fc ) = 
cos{(f)) — 

Afc = 
where 

A() = A© + 7T 



cos A 
cos(/3 fc ) 



2 < n v v < 
otherwise 



Mir 
2 



(16) 



(17) 



where A© is the position of the sun at the time the first 
ring, k is the ring number, 8 is the angular separation be- 
tween subsequent spin axis positions, n StP is the number of 
periods within 2n, and A is the amplitude. The values of 
these parameters used here are 8=2.5', n a =2, n p = 2.05, 
and yl=10° . This value of 8, given a repointing once per 
hour, keeps the spin axis directed away from the sun. 



5.2 The input point source catalogue 

The approximate numbers of point sources visible with 
Planck may be predicted by using th e IRAS point source 
catalogue (PSC, Beichm an et alJ <ll988f) ). This is a catalogue 
of some 250,000 well- confirmed point sources, providing po- 
sitions, flux densities at 12, 25, 60 and lOOixm, uncertain- 
ties and various cautionary flags which are given for each 
source. The selection of objects from this catalogue and the 
extrapola tion of their fluxes to Planck frequencies, is dis- 
cussed in lHarrison et all i2004t) . Also included in the input 
catalogue is the Wilkinson Mi crowave Anisotropy Probe, 
WMAP, point source catalogue. iBennett et al.l ([2003) . This 
consists of 208 extragalactic sources as seen in the WMAP 
maps, at Planck LFI frequencies. These sources may be ex- 
trapolated to Planck HFI frequencies using the spectral in- 
dices provided in the WMAP catalogue. The input point 
source catalogue constructed this way contains 5796 extra- 
galactic point sources and 2286 galactic sources. 



Table 1 . The number of galactic and extragalactic sources with 
fluxes greater than 30 times the nominal level of the point source 
sensitivity in the ring data. Numbers in italics correspond to a 
polarised detector pair. 

Freq No. with SNR ^ 30 

(GHz) Extragalactic Galactic Total 



30 




19 









19 


44 




9 









9 


70 




2 









2 


100 




17 




1 




18 


143 


24 


13 


2 


2 


26 


15 


217 


19 


m 


22 


u 


41 


n 


353 


13 


5 


99 


68 


112 


73 


545 


21 




258 




279 




857 


81 




654 




735 




Any 




110 


654 




764 





Table 2. The total intensity of a point source in the ring data, 
required for a 1 a detection, assuming the goal noise levels are 
attained. 



Freq 
(GHz) 


Iff noise 
Unpolarised detector 


(ring) (mjy) 

Polarised detector 


pair 


30 






147 


11 






230 


70 






346 


100 






102 


143 


85 




120 


217 


98 




139 


353 


184 




261 


545 


290 






857 


332 







Table Q shows the number of galactic and extragalactic 
sources detectable above a signal-to-noise ratio of 30 in the 
rings, for each of the Planck frequencies, under the assump- 
tion of the goal noise levels, shown in Table El 



5.3 Generating simulated data 

Instead of simulating time ordered data, TOD, we simulate 
directly the list of detections of source transits, as delivered 
by the Level2 DPC from their analysis of the TOD. This list 
of detections includes the position in phase and the ampli- 
tude observed for the transit, together with their respective 
errors. Additionally, the list includes the number of the de- 
tector which made the observation and the ring number on 
which it occurred. 

When a source will be observed by a detector depends 
upon when the line-of-sight of the detector passes close to 
the position of that source. This in turn depends upon the 
scanning strategy employed, the focal-plane layout and the 
values of the geometric calibrations parameters as discussed 
above. 

The simulations assume that the spin axis is reposi- 
tioned every hour and that the nutation effects are small so 
that the individual scans may indeed be combined to reach 
the nominal sensitivity to point sources in rings, as shown 
m Table H Any scanning strategy which meets this proviso 
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may be employed. The simulations also allow the variation 
in time of the values of the geometric-calibration parame- 
ters, as may be expected from the discussion in Section f3. 21 
As any time variation in the parameters is expected to 
be slow, as also discussed in Section 13.21 the parameters 
will be constant over the time frame of an individual ring. 
Hence, the instantaneous offset in a geometric-calibration 
parameter, 7, may be defined on each ring by: 



7 (r,) = 70 + 71 



(18) 



where Fi is the current ring number, T max is the final ring 
of the mission and the ring numbers start from zero. The 
systematic offset in the parameter, 70, is hence defined as 
the instantaneous offset of the parameter exactly half-way 
though the mission and the drift in the value of the pa- 
rameter, 71, is defined as the total drift in the value of the 
parameter over the course of the mission. 

For every spin axis position, the instantaneous values of 
the geometric-calibration parameters are assessed and used 
to reconstruct the lines-of-sight of the focal plane. If a source 
is located nearby, the amplitude with which the source is 
observed is assessed and if above a threshold signal-to-noise 
ratio it may be included in the simulated list of detections. 
Errors on the amplitude of the detections are generated as- 
suming the goal values of the noise in the ring and white 
noise. In order to generate realistic errors in the phase of 
the detection, an assessment must be made on how well the 
position of the transit in the scan (phase) direction may 
be measured. If the beams are Gaussian the position of the 
peak of the transit, ipt, may be found by minimising: 



Ai 



(19) 



where a - ;, is the beam sigma and Ai and ipi are the am- 
plitudes and phases of the i th sample in the transit. By 
simulating point source transits it is possible to establish 
a relationship between the peak amplitude for the transit 
and the error in the recovered phase for the transit. Hence, 
an empirical relationship may be found between the signal- 
to-noise ratio of the detection of the point source and the 
magnitude of the error in the recovered phase of the transit, 



(20) 



o> _ 1.7 
~ b ~ SNR 

where SNR is the signal-to-noise ratio of the detection. This 
expression was used to determine the magnitude of the error 
in the phase to be included in the simulations. This phase 
error will dominate any errors in the phase which result from 
those geometric-calibration parameters determined soley by 
the star tracker, such as errors in the velocity-phase relation. 



6 RESULTS 

Once a simulated list of detections has been generated, it 
may be analysed using the methods discussed in Section [I] 
The range of offsets in the geometric-calibration parameters 
which may be successfully recovered and the accuracies to 
which the recovered values may be attained may then be 
investigated. Unless otherwise stated, the simulated list of 



-100 -50 50 100 

input Si* REFo (orcsec) 

Figure 7. This figure shows the ability of this method to recover 
systematic offsets in the reference phase, 54> r ef 1 in the case of the 
sinusoidal scanning strategy. The differences between the input 
and recovered values of the systematic offset in the reference phase 
are plotted against the input values. For each input value, one 
hundred noise realisations were performed and the mean of the 
difference plotted with its error. The grey line is the global mean 
of these differences and the dashed lines enclose the region within 
\a of this global mean. The recovery of Sip re f can hence be seen 
to be unbiased. 



detections is generated using the goal noise levels in Ta- 
bleland only includes detections from extragalactic point 
sources above a threshold signal-to-noise ratio of 30. 

Recovering a pointing offset depends upon being able to 
determine which detections are due to an individual point 
source. The largest magnitude of pointing offset which may 
be recovered is therefore related to the beam sizes of the 
30 GHz channel, which has the largest beams. It has proved 
possible to successfully recover a pointing offset of ~ 20', 
which is greatly in excess of the expected magnitude of the 
offset, as discussed in Section [3.21 

Figure Q shows the difference between the value of the 
systematic offset in the reference phase, 5ip rB j , input to 
the simulations and the recovered value of 8ip re f against 
the input value, for the sinusoidal scanning strategy. This 
figure shows that the methods presented here have suc- 
cessfully recovered Stp re f over a range of values which ex- 
ceed those which may be expected. One hundred noise re- 
alisations per value of 5ip r ef investigated were performed, 
with each mean and its associated error plotted. The grey 
line shows the mean value over all values of <5VVe/ with 
the dashed lines representing the la errors in this global 
mean. Similar analyses have been performed for all the other 
geometric-calibration parameters for both scanning strate- 
gies. It is found that they may be successfully recovered over 
the range of possible offset values. The reference phase and 
roll angle are found to have an unbiased recovery, but the 
boresight angle is particularly sensitive to the errors in the 
recovered point source positions and so is vulnerable to small 
biases in the recovery of its parameters. These biases were 
found to be of the order of a few hundredths of an arcsecond 
to a tenth of an arcsecond depending on the scanning strat- 
egy and the signal-to-noise threshold used for inclusion of 
detections in the analysis. Using only the highest signal-to- 
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-100 



-50 50 

input 6a FRPo (arcsec) 



100 



Figure 8. This figure shows the difference between the input 
and recovered values of the systematic offset in the boresight an- 
gle, 8apnp Q , versus the values input to the simulations, which 
used the sinusoidal scanning strategy (top panel) or the preces- 
sional scanning strategy (bottom panel). For each input value, 
one hundred noise realisations were performed and the mean of 
the difference plotted with its error. The grey line is the global 
mean of these differences and the dashed lines enclose the region 
within la of this global mean. The recovery of SctpRp is seen to 
be biased at the level of a few hundredths of an arcsecond in the 
case of the sinusoidal scanning strategy and unbiased in the case 
of the preccssional scanning strategy. 



noise detections increases the likelihood of a biased recovery 
of these parameters as their assessment is then limited to 
a very few point sources. Given the pointing requirements, 
discussed in Section 13.31 a potential bias in the recovery 
of a geometric calibration parameter at this level is not a 
significant cause for concern. Figure |H] shows the recovery 
of the systematic offset in the boresight angle, SaFRp , for 
the sinusoidal and precessional scanning strategies, in the 
upper and lower panels respectively. In the case of the sinu- 
soidal scanning strategy the recovery of Scifrp was found 
to be biased at the order of a few hundredths of an arcsec- 
ond, whereas no bias in the recovered value was found when 
the experiment was repeated using the precessional scanning 
strategy. 

The simulations may also be used as a check on whether 
the errors found for the recovered values of the geometric- 
calibration parameters are a true representation of the un- 
derlying errors in each of these parameters. The simulations 
reveal that the dispersion in the mean recovered value of a 
parameter is consistent with the calculated error in the pa- 
rameter over the range of input values investigated, again 
with the exception of the boresight parameters in which the 
calculated error may be overestimated relative to the disper- 
sion in the recovered value. The size of this overestimation is 
also found to have some dependence on the scanning strat- 
egy employed and the threshold signal-to-noise ratio used. 

Figure |U] shows the dispersion in the mean recovered 
value of Soifrpo and the mean calculated error in the same 
against the value input to the simulations, which used the 
precessional scanning strategy. This shows that the calcu- 
lated error is representative of the actual error in the recov- 
ered value of the parameter. The error in the recovered offset 



-100 -50 50 100 

input <5a FRPo (arcsec) 

Figure 9. This figure shows the dispersion in the mean recovered 
value of SapRp , black crosses, and the mean calculated error in 
the value of &appp , grey crosses, against the value of SapRp Q 
input to the simulations, which used the precessional scanning 
strategy. The dashed lines enclosed the lcr region about the mean 
dispersion, and the grey line is the mean value of the calculated 
error. The calculated errors for 8ctppp are seen to be represen- 
tative of the underlying error in the recovered value of 5appp ■ 
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-200 -100 100 200 

input <5a FR p, (arcsec) 

Figure 10. This figure shows the dispersion in the mean recov- 
ered value of 8appp 1 , black crosses, and the mean calculated 
error in the value of SctpRp 1 , grey crosses, against the value of 
Scippp 1 input to the simulations, which used the precessional 
scanning strategy. The dashed lines enclosed the lcr region about 
the mean dispersion, and the grey line is the mean value of the 
calculated error. The calculated errors for 8appp 1 are seen to 
be overestimates of the actual error in the recovered values of 
SctFRP 1 , again there is no dependence of the errors on the value 
of 5a F Rp 1 



also has no discernible relationship with the input value of 
the offset. Figure 117)1 shows the overestimation of the error 
in the drift of the boresight angle, 8otpRp 1 , as compared 
to the dispersion in the recovered value of this parameter. 
Again the precessional scanning strategy is used and no de- 
pendence of the errors on the initial offset in the parameter 
is found. 
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Table 3. Comparing the errors in the recovered values of the 
geometric-calibration parameters in the cases of the sinusoidal 
and precessional scanning strategies. In both cases only extra- 
galactic point source detections with a signal-to-noise ratios of 30 
or greater were used in the analysis. 



Scanning Strategy: 


sinusoidal 


precessional 




(arcsec) 


(arcsec) 




0.22 


0.16 


°V>re/l 


0.17 


0.19 




0.20 


0.15 


eT °>FRP 1 


0.92 


0.73 


°"P0 


4.2 


3.6 




12.6 


14.5 



Table 4. This table shows the errors in the recovered values of 
the geometric-calibration parameters for different signal-to-noise 
ratio thresholds, for the inclusion of detections in the analysis, in 
the case of the sinusoidal scanning strategy. 



Threshold SNR: 


40 


30 


20 




(arcsec) 


(arcsec) 


(arcsec) 




0.25 


0.22 


0.21 




0.19 


0.17 


0.16 




0.24 


0.20 


0.17 


a c*FRP 1 


1.10 


0.92 


0.87 


°>a 


4.7 


4.2 


3.8 




14.3 


12.6 


11.3 



The disagreement of the errors and the occasional bias 
in the recovered values found for the boresight parameters 
are both due to errors in the recovered positions of the point 
sources used in this analysis. Simulations in which the cor- 
rect positions of the point sources are used show no bias in 
the recovered values and no inconsistencies between the cal- 
culated errors and the dispersions in the recovered values. 
This also explains the dependence on the scanning strategy 
used, as this affects the errors in the recovered positions of 
the point sources. Limiting the evaluation of the boresight 
parameters to the top four frequency channels, which have 
the smallest beams, is found to minimise the biases in the 
recovered values. If required, an assessment of whether the 
recovered values of the boresight parameters are likely to be 
biased may be made using the actual point sources observed 
and scanning employed by Planck. 

Table El compares the errors in the recovered values of 
the geometric calibrations parameters, when the sinusoidal 
and precessional scanning strategies are used. Due to the 
geometry of the ring crossings the positions of the point 
sources are attained to a slightly higher accuracy when the 
precessional scanning strategy is used, and this results in the 
slightly lower errors in the case of the precessional scanning 
strategy, especially in the case of the boresight parameters. 

Table |l]shows the calculated errors in the recovered val- 
ues of the geometric calibrations parameters for threshold 
signal-to- noise ratios of 40, 30 and 20. This shows that in- 
cluding the lower signal-to-noise ratio detections has very 
little impact on the accuracies to which the geometric- 
calibration parameters may be attained, in the case of 
ideal conditions. The simulations, however, may also be 



Table 5. The errors found in the recovered values of the 
geometric-calibration parameters using the sinusoidal scannning 
strategy and extragalactic source detections with signal-to-noise 
ratios greater than 30, where (i) uses the goal noise levels, (ii) 
excludes the 857 GHz frequency channel, (iii) doubles error in the 
phase of the peak of the transits and (iv) uses double the goal 
noise levels. 





(i) 


(") 


(iii) 


(iv) 




(arcsec) 


(arcsec) 


(arcsec) 


(arcsec) 




0.22 


0.25 


0.45 


0.64 




0.17 


0.27 


0.36 


0.45 




0.20 


0.36 


0.36 


0.69 


&a FRPl 


0.92 


1.56 


1.71 


2.96 


°>0 


4.2 


4.3 


8.9 


10.8 


°>1 


12.6 


12.9 


26.2 


36.1 



Table 6. The maximum pointing error found for an HFI and 
LFI detector using the errors in the recovery of the geometric- 
calibration parameters shown in Table 151 where (i) uses the goal 
noise levels, (ii) excludes the 857 GHz frequency channel, (iii) dou- 
bles error in the phase of the peak of the transits and (iv) uses 
double the goal noise levels. 





(i) 


(ii) 


(iii) 


(iv) 




(arcsec) 


(arcsec) 


(arcsec) 


(arcsec) 


HFI (545 GHz) 


0.62 


0.95 


1.20 


1.93 


LFI (44 GHz) 


0.86 


1.13 


1.74 


2.55 



used to investigate the errors in the recovered values of 
the geometric-calibration parameters, under non-ideal con- 
ditions. Table 03 shows the errors attained in the geometric- 
calibration parameters using all the detections with a signal- 
to-noise ratio above the threshold value of 30, for four differ- 
ent scenarios. The errors attained using the goal noise levels, 
shown in Table are compared against those using double 
these goal noise levels, as well as those where the error in 
the phase of the peak of each transit was doubled. Given 
that the majority of the point sources detectable by Planck 
appear only in the top frequency channel, the above analysis 
was performed without this channel, to investigate whether 
its loss would destroy our ability to acheive the pointing 
requirements. The errors in the geometric-calibration pa- 
rameters shown in Table |^| may be expressed as errors in 
the pointing reconstruction. The largest errors in the recon- 
structed pointing will occur at the beginning and end of the 
mission when the errors in the drift parameters will make 
their largest contributions, as may be seen from eauation llSI 
Table H3 shows the errors in the pointing reconstruction, as 
found for the errors in the geometric-calibration parameters 
shown in Table^J for the detectors in the HFI and LFI which 
have the largest uncertainty in their positions. This is due 
to their location in the focal plane with respect to the FRP, 
hence their greater sensitivity to errors in the roll angle. Ta- 
ble El shows that even in the case of doubling the goal noise 
levels, the pointing requirements are still easily achievable. 
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7 DISCUSSION 

The methods presented here show that recovering the 
geometric-calibration parameters as part of the initial stages 
of the construction of the final source catalogue, successfully 
extracts any offsets in the values of these parameters, using 
only the detections due to extragalactic point sources. The 
accuracy to which the geometric-calibration parameters may 
be attained far exceeds the pointing accuracy requirements 
as discussed in sect ion |3~31 so much so that it has proved pos- 
sible to recover the geometric-calibration parameters under 
the non-ideal conditions of the loss of a frequency channel, 
and double the goal noise levels. The accuracies to which the 
geometric-calibration parameters may be recovered, are such 
that the errors in the pointing reconstruction due to errors 
in the geometric-calibration parameters are of the same level 
as the errors in the pointing reconstruction due to the uncer- 
tainties in the mean spin axis position recovered by the star 
tracker which are expected to be of the order of l"-2". The 
achievable accuracies for the geometric-calibration parame- 
ters as found above did not include errors in the remaining 
geometric-calibration parameters found solely from the star 
tracker such as errors in the mean spin axis position or the 
velocity-phase relation. The expected level of these errors, 
however, has a negligible effect on the errors found above. 

The focal-plane layout, relative to the FRP, was as- 
sumed to be fixed and known. If thermal variations produce 
offsets in the positions of the detectors relative to the FRP, 
these will also need to be recovered using the science data. 
It is anticipated that the any positional offsets in the focal- 
plane layout will be recovered using the planetary transits 
of the focal plane. 

The methods presented here have assumed Gaussian 
beams. In reality, however, the Planck beams will not be 
Gaussian. If the beams are not symmetric about the scan 
direction then the phase, corresponding to the peak ampli- 
tude of the transit, will depend on the ordinate of the point 
source. Any dependence of the phase, of the transit, on the 
ordinate of the point source must be included in the analysis, 
and may lead to increases in the errors of the phases for the 
detections. We have however demonstrated that doubling 
the error in the phase of each transit has the effect of roughly 
doubling the error in the pointing reconstruction due to the 
geometric calibrations parameters, and these resultant er- 
rors still easily meet the required pointing accuracy. There 
is therefore plenty of scope for increased uncertainties in the 
positions of the transits. 

Due to the number of available bright point sources, the 
time resolution of the recovery of the geometric-calibration 
parameters is poor and any fast evolution in a parameter 
will not be recoverable from the science data. Of concern is 
whether an unsolved variation in a parameter could bias the 
recovery of the parameter through an uneven distribution 
of detections. To investigate this possibility, the effect of an 
unsolved drift in the parameters was investigated. An un- 
recovered drift was found not to affect the recovered value 
of the systematic offset until the pointing requirements are 
exceeded by the unrecovered drift itself. It is therefore likely 
that any unsolved for variations will not affect the recov- 
ery of the systematic offsets and drifts while they are small 
enough not to be a problem in and of themselves. 
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